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Abstract. In a recent paper [Yj we give the first rigorous derivation of the cel- 
ebrated Ginzburg-Landau (GL) theory, starting from the microscopic Bardeen- 
Cooper-Schrieffer (BCS) model. Here we present our results in the simplified 
case of a one-dimensional system of particles interacting via a <5-potcntial. 



1. Introduction and Main Results 

1.1. Introduction. In 1950 Ginzburg and Landau [T] presented the first satis- 
factory mathematical description of the phenomenon of superconductivity. Their 
model examined the macroscopic properties of a superconductor in a phenomeno- 
logical way, without explaining its microscopic mechanism. In the GL theory the 
superconducting state is represented by a complex order parameter ip(x), which 
is zero in the normal state and non-zero in the superconducting state. The order 
parameter ip( x ) can he considered as a macroscopic wave-function whose square 
|V>(aO| 2 is proportional to the density of superconducting particles. 

In 1957 Bardeen, Cooper and Schrieffer 2 formulated the first microscopic ex- 
planation of superconductivity starting from a first principle Hamiltonian. In a 
major breakthrough they realized that this phenomenon can be described by the 
pairing-mechanism. The superconducting state forms due to an instability of the 
normal state in the presence of an attraction between the particles. In the case of 
a metal the attraction is made possible by an interaction through the lattice. For 
other systems, like superfluid cold gases, the interaction is of local type. In the BCS 
theory the superconducting state, which is made up by pairs of particles of opposite 
spin, the Cooper-pairs, is described by a two-particle wave-function a(x,y). 

A connection between the two approaches, the phenomenological GL theory and 
the microscopic BCS theory, was made by Gorkov [3] who showed that, close to the 
critical temperature, the order parameter ip(x) and the pair-wavefunction a{x,y) 
are proportional. A simpler argument was later given by de Gennes [4] . 

Recently we presented in a mathematical proof of the equivalence of the 
two models, GL and BCS, in the limit when the temperature T is close to the 

1/2 

critical temperature T c , i.e., when h = {(T c — T)/T c ] 1 <C 1, where T c is the critical 
temperature for the translation-invariant BCS equation. The mathematical aspects 
of this equation where studied in detail in [HI El El [101 E] ■ In the present paper 
we present this result in the simplified case a of one-dimensional system where the 
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particles interact via an attractive contact interaction potential of the form 

V(x - y) = -aS(x — y) witha>0. (1.1) 

We assume that the system is subject to a weak external potential W, which 
varies on a large scale 1/h compared to the microscopic scale of order 1. Since 
variations of the system on the macroscopic scale cause a change in energy of the 
order h 2 , we assume that the external potential W is also of the order h 2 . Hence 
we write it as h 2 W(hx), with x being the microscopic variable. The parameter h 
will play the role of a semiclassical parameter. 

We will prove that, to leading order in h, the Cooper pair wave function a(x, y) 
and the GL function ip(x) are related by 

a(x,y)=^(h^p^a°(x-y) (1.2) 

where a is the translation invariant minimizer of the BCS functional. In particular, 
the argument x of the order parameter ip(x) describes the center-of-mass motion 
of the BCS state, which varies on the macroscopic scale. To be precise, we shall 
prove that a(x, y) — ^(tp(hx) + ip(hy))ct° {x — y) to leading order in h, which agrees 
with p. 21) to this order. 

For simplicity we restrict our attention to contact potentials of the form (|1.1[) . 
but our method can be generalized to other kinds of interactions; see [7] for details. 
The proof presented here is simpler than the general proof in |7| which applies 
to any dimension d < 3. There are several reasons for this. First, there is no 
magnetic field in one dimension. Second, for a contact interaction the translation 
invariant problem is particularly simple and the corresponding gap equation has an 
explicit solution. Finally, several estimates are simpler in one dimension due the 
boundedness of the Green's function for the Laplacian. 

1.2. The BCS Functional. We consider a macroscopic sample of a fermionic 
system, in one spatial dimension. Let fi G M denote the chemical potential and 
T > the temperature of the sample. The fermions interact through the attractive 
two-body potential given in (II. ip . In addition, they are subject to an external force, 
represented by a potential W{x). 

In BCS theory the state of the system can be conveniently described in terms of 
a 2 x 2 operator valued matrix 




satisfying < T < 1 as an operator on L 2 (R) L 2 (R). The bar denotes complex 
conjugation, i.e., a has the integral kernel a(x, y). In particular, T is assumed to 
be hcrmitian, which implies that 7 is hermitian and a is symmetric (i.e, 7(2, y) — 
"/(y,x) and a(x,y) = a(y,x).) There are no spin variables in T. The full, spin 
dependent Cooper pair wave function is the product of a with an antisymmetric 
spin singlet. 

We are interested in the effect of weak and slowly varying external fields, de- 
scribed by a potential h 2 W(hx). In order to avoid having to introduce boundary 
conditions, we assume that the system is infinite and periodic with period h^ 1 . In 
particular, W should be periodic. We also assume that the state T is periodic. The 
aim then is to calculate the free energy per unit volume. 
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We find it convenient to do a rescaling and use macroscopic variables instead of 
the microscopic ones. In macroscopic variables, the BCS functional has the form 

J- BCS (r) :=Tr(-h 2 \7 2 - ii + h 2 W(x))j-TS{T)-ah J \a(x,x)\ 2 dx (1.3) 

where C denotes the unit interval [0,1]. The entropy equals S(T) — — TrTlnT. 
The BCS state of the system is a minimizer of this functional over all admissible T. 

The symbol Tr in (|1.3p stands for the trace per unit volume. More precisely, if B 
is a periodic operator (meaning that it commutes with translation by 1), then TrB 
equals, by definition, the (usual) trace of X-B, with \ the characteristic function 
of C. The location of the interval is obviously of no importance. It is not difficult 
to see that the trace per unit volume has the usual properties like cyclicity, and 
standard inequalities like Holder's inequality hold. This is discussed in more detail 
in[ZJ. 

Assumption 1. We assume that W is a bounded, periodic function with period 1 
and J c W(x) dx = 0. 

1.2.1. The Translation-Invariant Case. In the translation invariant case W = one 
can restrict J rBCS to translation invariant states. We write a general translation 
invariant state in form of the 2x2 matrix 

/ 7HW) a (-ifeV) \ 

v &{-ihV) l-T(-iW) J ' [ ' 

that is, 7 = [r]n and a = \T]i2 have integral kernels 

7 (ar, y ) = 2- I j(h P y p{x - y) dp and a(x, y) = ^- f &(hp)e tp{x - y ^ dp . 

27T ./» 27T ./» 



The fact that T is admissible means that a(p) = a(—p), that < j(p) < 1 and 
|5(p)| 2 < 7(p)(l ~ t( — p)) f° r an y P S R- For states of this form the BCS functional 
becomes 



T^D= I (h 2 p 2 -^(hp)^-~T [ S (t(hp))^-ah 



a(hp)^- 

Z7T 



2 



(1.5) 



with S{T (p)) = — Trc2 t(p) lnf (p) and f (p) the 2x2 matrix obtained by replacing 
-iV by p in (fT3]) . 

In the following, we are going to summarize some well-known facts about the 
translation invariant functional (|1.5[) . For given a > 0, we define the critical tem- 
perature T c > by the equation 

i = n^M^w^dp 

a Jr P 2 - H 2?r ' 

The fact that there is a unique solution to this equation follows from the strict 
monotonicity of t/tanht for t > 0. If T > T c , then the minimizer of (|1.5[) satisfies 
q = and j(hp) = (1 + exp((/i 2 p 2 - ^/T))- 1 . If < T < T c , on the other hand, 
then there is a unique solution Ao > of the BCS gap equation 

1 f 1 dp 



(1.7) 
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where 

tanh ^/(p2 - M )2 + Ag) 



HIM 



Moreover, the minimizer of (jl.5|l is given by 

f°(/ip) = (l + exp(i J ffi o (M)) _1 (1-9) 

with 

P 2 - M -A 
-A -p 2 + /u 

Writing a (ftp) = [r°(/ip)]i2 one easily deduces from (jl.9|) that 

*°w=4fer (L10 » 

To summarize, in the case TU = the functional J-" BCS has a minimizer T° for 
< T < T c whose off-diagonal element does not vanish and has the integral kernel 

We emphasize that the function a depends on T. For T close to T c , which is the 
case of interest, we have Ao ~ const (1 — T/Tc) 1 ^ 2 . 

1.3. The GL Functional. Let t/j be a periodic function in H^ oc (M.). For numbers 
&i, 63 > and &2 G K the Ginzburg-Landau (GL) functional is given by 

£ty) = jT (&i|^'(x)| 2 + 6 2 ^(x)|V(a;)| 2 + 6 3 (l - |^(a;)| 2 ) 2 ) dx . (1.12) 

We denote its ground state energy by 

£ GL =inf{£ty)|^ 

Under our assumptions on W it is not difficult to show that there is a corresponding 
minimizer, which satisfies a second order differential equation known as the GL 
equation. 

1.4. Main Results. Recall the definition of the BCS functional J" BCS in ([T3]). We 

define the energy F BCS (T, /z) as the difference between the infimum of F BCS over 
all admissible T and the free energy of the normal state 



(1.13) 



7o 
I-70 

with 70 = (1 + e (-k 2 v 2 +fc 2 w(*)-^)/T)-i. That iSj 

F BCS (T; m) = ^ ^BCS (r) _ ^-BCS ( p o) _ (L 14) 

Note that 

J- BCS (r ) = -T Trln (l + exp (- (-ft 2 V 2 - n + h 2 W \x))) /T) . (1.15) 

For small h this behaves like an (explicit) constant times h~ l . Under further 
regularity assumptions on W, (II . 15[) can be expanded in powers of h. We do not 
need this, however, since we are only interested in the difference F BCS (T, jj). 
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Since Tq is an admissible state, one always has F BGS (T,[i) < 0. If the strict 
inequality F BCS (T, /i) < holds, then the system is said to be in a superconducting 
(or superfluid, depending on the physical interpretation) state. 

THEOREM 1. Let Assumption [7] be satisfied, and let T c > be the critical 
temperature in the translation invariant case, defined in fl.6\) . Let D > 0. Then 
there are coefficients b\, b 2 and 63, given explicitly in (|1.20|) - (|1.22j) below, such that 

F BCS (T C (1 - Dh 2 ), p) = h 3 (E GL - 63) + o(h 3 ) (1.16) 

as h — > 0. More precisely, the error term o(h 3 ) satisfies 

- const h 3+ 3 < o(h 3 ) < const h 5 . 

Moreover, ifT is an approximate minimizer of J- BCS atT = T c (l — h 2 D), in the 
sense that 

j- BCS (r) < j- BCS (r ) + h 3 (E GL - 63 + e) (1.17) 

for some small e > 0, then the corresponding a can be decomposed as 

a(x, y) = ±(ip(x) + t/>(y))a Q (h- 1 (x - y)) + a(x, y) (1.18) 

with £ GL (ifj) < E GL + e + const h\ , a a defined in (TU\) . and 



\a(x,y)\ 2 dxdy < const ft, 3+ 3 . (1.19) 



Cxi 



1.5. The coefficients in the GL functional. In order to give explicit expressions 
for the coefficients in the GL functional we introduce the functions 

, . tanh(z/2) , . e 2z -2ze z -l , , 2e z (e z - 1) 

9o z = , 9i{z) = — , g 2 {z) = — -3-. 

z z z (l + e z Y z(e z + l) 

Setting, as usual, (3 C = T^ 1 we define 

e _ 2jg [g ((3 c (q 2 - fx)) - (3 c (g 2 - M )gi(/3c(g 2 - //))] dq 
or gi(Pc{g 2 -n)) , 

The three coefficients of the G-L functional turn out to be as follows, 

h = cD ^ Jj 9l (l3 c ( q 2 - f i)) + 2f3 c q 2 g 2 (M<l 2 - V))) ^, (1-20) 

b 2 = cD^- JjMq 2 -^)^. (1.21) 

and 

h = {cDf Pl ( gx(W-M)) ^ (L22) 

We shall now discuss the signs of these coefficients. First note that go(z) — 
zg\(z) = (zgo(z))' > and g\(z)/ z > 0, which implies that c > 0. Using g\[z)j z > 
again, we see that &3 > 0. In contrast, the coefficient b 2 may have either sign, 
depending on the value of /3 c /x (which depends on a and //). The coefficient fei is 
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again positive, as the following computation shows: using the fact that 32(2) 
g[(z) + (2/z)gi(z) we find 

7 n 01 f ( (o r 2 \n , fl 2 / //o / 2 xx , 2 5 i(/? c (g 2 -//))^ d 9 
&i=cD— / 5i(/3 c (g -~n)) + 2p c q 3i(p c (g - w) + 



2 9i(Pc{q 2 - (J>))\ dq 



cD^J^ (W(<z 2 - m)) + (gi(/3 c (g 2 - m))) + V 



</ 2 — /i / 27T 



cD !% [ q 2 gi(/3c(g 2 ~M)) ^ 



The claimed positivity is now again a consequence of g\(z)/z > 0. 



2. Sketch of the proof 



In the following we will consider temperatures T = T c (l — Dh 2 ). It is not difficult 
to see that the solution A of the BCS gap equation (|1.7p is of order A = 0(h). 

It is useful to rewrite the BCS functional in a more convenient way. Define A to 
be the multiplication operator 

A = A(x) = -i/)(x)Ao, 

where A is the solution of the BCS equation (| 1 . T[) for temperature T, and ip a 
periodic function in H 2 oc (M.). Define further 



-h 2 V 2 



fi + h 2 W(x) 

A h 2 V 2 



A 

fi - h 2 W(x) 



(2.1) 



Formally, we can write the BCS functional as 

1 



T U ^(T) 



Tr(-/i 2 V 2 
1 



(j, + h A W) + - TvH A T - TS(T) 



\tp(x)\ 2 dx — ha 



A iP(x) 



2ha 



a(x, x) 



dx. (2.2) 



The first two terms on the right are infinite, of course, only their sum is well-defined. 
For an upper bound, we can drop the very last term. The terms on the first line 
are minimized for Ta = 1/(1 + C A ), which we choose as a trial state. Then 



f bcs (t,v) < J- BCS (r A ) - j- BCS (r ) 

T 



(2.3) 



< 



■Tr 



ln(l 



e~T HA ) -kx(l + e~T Ho ) 



Aha 



\ip(x)\ 2 dx. 



To complete the upper bound, we have to evaluate Tr[ln(l + e~ ffA//T ) — ln(l + 
e~ H °/ T )]. This is done via a contour integral representation and semiclassical types 
of estimates. 

The lower bound is divided into several steps. We first aim at an a priori bound 
on a for a general state T, which has lower energy than the translation-invariant 
state. With 

o _ / -h 2 V 2 -fi A 
H ^>-\ A -/i 2 V 2 + /z 
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we can rewrite the BCS functional in the form 



j BU5 (r) 



Tr(-/i 2 V 2 -fi + h 2 W) + - tr H A( Y ~ TS(T) 



h 2 TrWj+-^-AZ-ha [ 
Aha J c 

d the definit: 
* (0)= 2^X 2A "v'W'.) 



2h-a- a(x ' X ' 



dx. (2.4) 



From the BCS equation and the definition of a in (jl.lip we conclude that 

Ao 



dp= 2aH> 



(2.5) 



and hence 

J- BCS (r)-^ BCS (r°) > ^H(r,r ) + / l 2 TrT^( 7 -7 )-a/i / \a(x,x) -a°(0)\ 2 dx, 

2 Jc 

(2.6) 

where % denotes the relative entropy 

H(r,r°) = | Qtr#£ o r - ts(t) + \tr h° a y° - ts(y": 

= Tr [r (lnT - lnT ) + (1 - T) (ln(l - T) - ln(l - T ))] . (2.7) 

Note that the left side of (I2.6[) is necessarily non-positive for a minimizing state T. 

One of the essential steps in our proof, which is used on several occasions, is 
Lemma [TJ This Lemma presents a lower bound on the relative entropy of the form 

U(T,T°) > Yt\h° (F-F ) 2 

i (Trr(i-r)-Trr (i-r°)) 2 

3|TrF(l-F)-Trr°(l-F°)|+TrF (l-r )' ( ' ' 

where H° = (1 - 2r r 1 ln((l - F°)/F°). In our case here, it equals K^(~ihV)/T, 
with defined in (ll.8[) . From (|2.6I) we deduce that for a minimizer Y 

> J" BCS (r) - J- BCS (F°) > TrK°(l - 7 ) 2 + h 2 Tr W(j - 7 °) 

+ J (<*(-, y)- a\^-)\K° T {-ihV) - a5{^)\a{ ■ , y) - a°(^)> 

1 r(Tr[7(l-7)-7 (l-7 )-|«| 2 + l«°| 2 ]) 2 

3|Tr[7(l-7)-7 (l-7 )-|a| 2 + |a | 2 ]|+Tr[7 (l-7 )-|a°| 2 ]' 1 ' j 

where (-|-) denotes the inner product in L 2 (M). Observe that the term in the second 
line is a convenient way to write Yr Kj.{a — a°)*(a — a Q )~ ah J c \a(x, x) — a°(0)\ 2 dx. 
From the first line on the right side and the Schwarz inequality together with the 
fact that K^-a8> 0, we obtain first that Tr if£( 7 - 7 ) 2 < 0(h 3 ). Together with 
the last line this further gives the a priori bound Hall! < 0(h). 

Next, we use that — aS has a as unique zero energy ground state, with 
a gap of order one above zero, and we can further conclude from (|2.9[) that a is 
necessarily of the form 

a(x,y) = \{^(x) + ^{y))a°{{x-y)/h)+P{x,y), 

with 1 1/3 1 1 2 < 0(h 3 ). This information about the decomposition of a then allows us 
to deduce, again by means of a lower bound of the type (|2.8|) . that the difference 
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J rBCS (r) — J rBCS (rA) is very small compared to h 3 . This reduces the problem to 
the computation we already did in the upper bound. 

3. Semiclassics 

One of the key ingredients in both the proof of the upper and the lower bound 
are detailed semiclassical asymptotics for operators of the form 

( -W-v + h*W(x) A( X ) \ 

Ua -{ A (x) h 2 V 2 +^-h 2 W{x) ) ■ [ > 

Here A(x) = —hi]j{x) with a periodic function ip, which is of order one as h — > 
(but might nevertheless depend on h). We are interested in the regime h — ¥ 
0. In contrast to traditional semiclassical results [121 116] we work under minimal 
smoothness assumptions on ip and W. To be precise, we assume Assumption [1] for 
W and that ip is a periodic function in H 2 oc (M.). 
Our first result concerns the free energy. 

THEOREM 2. Let 

/(*) = - hi (1 + e-*) , (3.2) 

and define 

, , f'(-z) - f'(z) tanh(iz) , 
9o(z) = 1 ' J 1 ' = , (3.3) 



ffi(«) = -ffo(*) = 12 + : = .2 f1le ,« ( 3 - 4 ) 



/'(-*) - /'(«) , /"(-«) + f"(z) _ e 2z - 2ze 

z 2 {l + e z ) 
and 

n ,,,2 /"'(s) - f" (-z) 2e' (e z - 1) 
52(«) =5i(«) + -fl'iW = = — -73 ■ (3-5) 

z z z (e z + 1) 

Then, for any ft > 0, 

| Tr [/(/Lff A ) - /(£#„)] = ft 2 ^ + /j 4 £ 2 + 0(^ 6 ) (||^||^ (c) + H\\ 2 h H c)) , (3.6) 
where 

2 JR 27T 

and 

E2 =-t-W'\\t I (9i(P(<l 2 A*)) + 2/? 2 <? 2 (/3( 2 - M))) 



2tt 

+ y(^IW>^3i(/% 2 ~M))§ 

More precisely, we claim that the diagonal entries of the 2x2 matrix- valued 
operator f(fiH&) — f((3H ) are locally trace class and that the sum of their traces 
per unit volume is given by (|3.6I) . We sketch the proof of Theorem [2] in Subsection 
6.21 below and refer to [7] for some technicalities. 

Our second semiclassical result concerns the behavior of (1 + exp(/3i/A)) _1 in 
the limit h — >■ Q. More precisely, we are interested in [(1 + exp(/3_ffA)) _1 ]i2, where 
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[ • ]i2 denotes the upper off-diagonal entry of an operator- valued 2x2 matrix. For 
this purpose, we define the H 1 norm of a periodic operator n by 

\\n\\ 2 Hl =Tr[ V *(l-h 2 V 2 )n] . (3.7) 

In Subsection 16.31 we shall prove 



THEOREM 3. Let 

p(z) = (l + e z )- 1 (3-8) 

and let go be as in (|3.3[) . Then 

WH A )} 12 = ^ (i,(x) g ((3(-h 2 \7 2 - //)) + goM-tfV 2 - fx)) ^(x)) + m +n 2 
where 

hi\\m<Ch 5 U\\ (3.9) 

and 

\\V2\\^<Ch 5 [u\\ 2 HHC) + n\\m { c)) ■ (3-10) 
4. Upper Bound 

We assume that T = T c (l — Dh 2 ) with a fixed D > and denote by Ao the 
solution of the BCS gap equation (|1.7[) . In the following we write, as usual, j3 = 
T" 1 = f} c (\ - Dh 2 )- 1 with f} c = T~ l . It is well known that the Ginzburg-Landau 
functional has a minimizer -0, which is a periodic H 2 oc {R) function. We put 

A(x) = -A ^(x), 

and define H A by ([2TT]) . 

To obtain an upper bound for the energy we use the trial state 

r A = (i + e ^)- 1 . 

Denoting its off-diagonal element by a a = [Ta] 12, we have the upper bound 
^ BCS (r A ) - ^ BCS (r ) = -i-Tr [ln(l + e~^) - ln(l + e"^°)] 



CXA\X, X) 



2ha 



dx (4.1) 



< - JL Tr [ln(l + e-^) - ln(l + e"^)] + ^UWl ■ 

The first term on the right side was evaluated in Theorem[2j Applying this theorem 
with replaced by (Ao/h)ip we obtain that 

j- BCS (r A )-^ BCS (r ) 
<- k 4^U\\llM-,)) dq 



4 h?" Y " 2 J M »"^^ tn 2vr 

+ (4.2) 



10 



R.L. FRANK, C. HAINZL, R. SEIRINGER, AND J. P. SOLOVEJ 



In the estimate of the remainder we used that ip is H 2 and that Ao < Ch. 

Next, we use that by definition (|1.7j) of Ao the first and the last term on the 
right side of (|4.2[) cancel to leading order and that one has 



h P A 0il,.ii2 f ( „ ta [7~2 ,A2 , a 2 \ „ rarJZ ..\\\ dc l 



8 ft, 4 7 R g 2 ~M 2tt 

We conclude that 

^ BCS (r A )-^ BCS (r ) 



ft 3 
< — 

~ 2 



^M 2 /" 

4 ft^H 



ffi(/3(<7 2 -M)) rf 9 



0(/i 5 ). (4.3) 



— /i 27T 

Up to an error of the order 0(h 5 ) we can replace j3 = /3 C (1 — Dh 2 ) -1 by /3 C on the 
right side. Our last task is then to compute the asymptotics of Ag/h. To do so, we 
rewrite the BCS gap equation (|1.7p as 

/3 C / 9o((3 c (q 2 - M)) ^ = - = /? / 3o(/?a/(9 2 -m) 2 +A 2 ) ^ . 
A simple computation shows that 

Aq = Dh 2 k[aom<i 2 -^-Pjf : ^mq 2 ~m dg , a } 

PcJr 2(< ? 2 - ai ) "y 

Inserting this into (I4.3[) and using the fact that £(V0 = E Gh we arrive at the upper 
bound claimed in Theorem [1] 

5. Lower Bound 

5.1. The relative entropy. As a preliminary to our proof of the lower bound, 
we present a general estimate for the relative entropy. In this subsection H° and 
< r < 1 are arbitrary self-adjoint operators in a Hilbert space, not necessarily 
coming from BCS theory. Let L° := (l + cxp(/3iJ )) . It is well-known that 

H(T, r°) = Tr {pH a T + LlnL + (1 - L) ln(l — L) + In (1 + exp(-jftffo))) 

is non- negative and equals to zero if and only if L = L°. Solving this equation 
for H°, i.e., H° = /3 _1 (ln(l - L°) - lnL°), we can rewrite H(T,T Q ) as a relative 
entropy, 

H(T,T°) = Tr [r (lnL - lnL°) + (1 — L) (ln(l -T) — ln(l — L ))] . (5.1) 

The following lemma quantifies the positivity of % and improves an earlier result 
from [5]. 
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Lemma 1. For any < T < 1 and any T of the form T° = (1 + e 0H ) 
H(T,T°) > Tr 



tanh(^J?°/2) 

i (Trr(i-r)-Trr (i-r )) 2 
+ 3 |Trr(i - r) - Trr°(i - r°)| + Trr°(i - r°) ' 

Proof. It is tedious, but elementary, to show that for real numbers < x, y < 1, 

, x 1-x 3 1 (x(l - x) - y(l - y)f 

xm — h (1 — x)m- >- — (x — y) + 



V 1-y l-2y 3 \x(l — x) — y(l — y)\ + y(l — y) 

Using joint convexity we see that 

(x(l -x)- y(l - y)f 



\x{\ -x)-y(l-y)\ + y(l - y) 

= 4 sup [6(1 - b) \x(l - x) - y{\ - y)\ - b 2 y(l - y)] . 

0<b<l 

Let us replace on the right side the modulus \a\ by max{a, — a}, and then use 
Klein's inequality (TSJ Section 2.1.4] for either of the expressions. This implies the 
result. □ 

5.2. A priori estimates on a. We begin by briefly reviewing some facts about 
the translation-invariant case W = 0; see also Subsection 11.2.11 Recall that T° 
denotes the minimizer of T in the translation-invariant case. It can be written as 
r° = (1 + e 0H ° & o)- 1 with 

o _ ( -h 2 V 2 - /i -A 
Ua « ~ { -A /i 2 V 2 + fx 

Here A is the solution of the BCS gap-equation (fLT]) and = T = T c (l - Dh 2 ). 
Notice the distinction between T° and To which was defined in (|1.13|) . The latter 
one, To, contains the external potential W and has no off-diagonal term. 

Recall also that we denote the kernel of the off-diagonal entry a = [r°]i2 by 
a°((x — y)/h), which is explicitly given in (jl.lll) . From this explicit representation 
and the fact that Ao < Ch we conclude, in particular, that 

\\a°\\l= [ dy [ dx\a°((x-y)/h)\ 2 =h [ \a° (x)\ 2 dx < C h . (5.2) 



Moreover, the BCS gap-equation JT77J) is equivalent to 

(K^(-ihW) - ah6(x))a°(x/h) = . (5.3) 
This implies, in particular, that 

Tr K%(-ihV)a°aP - ah\a°(0)\ 2 = 0. (5.4) 

Now we turn to the case of general W. Our goal in this subsection is to prove 
that the a of any low-energy state satisfies bounds similar to (|5.2[) and ([57 



Proposition 1. Any admissible T with T BCS (T) < J" BCS (r°) satisfies 

\\a\\l = [ dx [ dy\a(x,y)\ 2 <Ch (5.5) 



n 
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<TrK%(-ihV)aa- ah J \a(x, x)\ 2 dx < Ch 3 , 



(5.6) 



where a = \T]i%. 

Proof. We divide the proof into two steps. 

Step 1. Our starting point is the representation 

T BGS {r)~T BCS {r ) = ±H(T,r°) + h 2 Tr~fW-ah f \a(x, x) - a°(0)| 2 dx 

(5.7) 

for any admissible T, with the relative entropy 7i(T,T°) defined in (|5.1[) . We note 
that T° is of the form (1 + e^") -1 with H° = H% q . We use Lemma[T]to bound 
T-L(T,T°) from below. Since x i— > x/tanhx is even, we can replace H\ by its 
absolute value E(—ihV) = (— h 2 V 2 — fi) 2 + A 2 ,, and thus 



Tr 



An 



tanh § H^ o 



(r-r ) 2 



Tr[K°(-ih\7)(T~r Q ) 2 ] 



with 



2Tr K%(-ihV)(i - 7 ) 2 + 2TriT£(-iW)(a - a ) (a - a ) 

£(M 



tanh 



PE{hp) 



from (TEH]). With the aid of Lemma [Hand the assumption J rBCS (r) < J" BCS (r°) 
we obtain from ()5.7|) the basic inequality 



>Tr^(-i/iV)( 7 -7 ) 2 + h 2 TrjW 

+ Tr K%.(-ihV)(a - a°){a - a ) - ah I \a(x, x) - a°(0)| 2 dx 



+ 



T 


7(1 - 7) - 7°(1 - 7 ) - aa + a°a° 


r 


3 


Tr 


7(1 — 7) — 7°(1 — 7 ) — aa + a°a° 


+ Tr 


7 °(1 - 7 ) - a°a° 

( 



(5.8) 

In the following step we shall derive the claimed a priori estimates on a from this 
inequality. 

Step 2. We begin by discussing the first line on the right side of (|5.8[) . Using the 
fact that Tr Wj° = (since W has mean value zero) and the Schwarz inequality 
we obtain the lower bound 

Tr lf£(-iW)(7 - 7°) 2 + h 2 TrW-/ 

= \ TrK° HW)( 7 - 7 ) 2 + \ TrX° (7 - 7 ) 2 + h 2 Tr W(~f - 7 ) 

> ± Tr J^-i/iV) (7 -7 ) 2 - ^iTrW (K^i-ihV))^ 1 W 



> iTr^(-iW)( 7 -7 ) 2 



2 

Ch 3 . 



(5.9) 



The last step used that Tr W (K%,(-ihV)) 1 W < \\ Tr K^-ihV)- 1 < Ch,- 1 . 
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Next, we treat the second line on the right side of (|5.8p . Recall that the BCS 
gap equation in the form ()5.3[) says that the operator Kj,(—ihV) — ahS(x) has an 
eigenvalue zero with eigenfunction a (x/h). Hence 



Tr K%(-ihV)(a - a°)(a - a ) - ah J \a{x, x) - a°(0)| 2 dx 
= Tr K^(-ihV)aa - ah J \a(x, x)\ 2 dx . 

Since a delta potential creates at most one bound state, zero must be the ground 
state energy of Kj,(— iftV) — ahS(x), and we deduce that 



Tr K^(-ihV)aa - ah J \a(x, x)\ 2 dx > . 



This information, together with (|5.9I) and (|5.8I) . yields 

Tr(l-/i 2 V 2 )( 7 -7°) 2 < CTr^(-iW)(7-7 ) 2 < Ch 3 



Tr K^(-ihV)aa -ah / \a(x, x)\ 2 dx < Ch 3 
Jc 



(5.10) 
(5.11) 



and 



(Tr 


7(1 - 7) - 7°(1 - 7 ) - aa + a°a° 


) 


2 


Tr 


7(1 — 7) — 7°(1 — 7°) — oloi + a°a° 


+ Tr 


7°(1- 


- 7 ) - 



< Ch 3 



_ (5.12) 
7°(1 — 7 ) — a°a° < Ch^ 1 from the explicit solution in the 



We know that Tr 
translation invariant case, and therefore (|5.12[) yields 



Tr 



7(1 - 7) - 7°(1 - 7 ) - aa + a°a° 



< Ch. 



(5.13) 



In order to derive from this an a priori estimate on a we use (|5.10l) and the Schwarz 
inequality to bound 



|Tr(7-7°)| < /r 2 Trfr£(-i/iV)(7-7 ) 2 + h 2 Tr (K%(-ihV)) 1 < Ch 



and 



|Tr( 7 -{rY)\ = |Tr( 7 -7 U )( 7 + 7 U )| < h~ J Tr( 7 - 7 U )^ + h* Tr( 7 + j v ) z < Ch . 

Finally, since Tra°a° < Ch (see (|5.2p ). we conclude from (|5.13[) that Traa < Ch, 
as claimed. □ 

5.3. Decomposition of a. Here we quantify in which sense a(x, y) is close to 
i (ip(x) + ip(y)) a°(h~ 1 (x — y)). There is one technical point that we would like to 
discuss before stating the result. The asymptotic form i (ip(x) + ip(y)) a°(h~ 1 (x — 
y)) will allow us in the next subsection to use the semiclassical results in a similar 
way as in the proof of the upper bound. Our semiclassics, however, require ip to be 
in H 2 . While we naturally get an H 1 condition, the H 2 condition is achieved by 
introducing an additional parameter e > 0, which will later chosen to go to zero as 
h^O. 
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Proposition 2. Let T be admissible with J rBCS (r) < J" BCS (r°). Then for every 
sufficiently small e > h > 0, the operator a — [T]\2 can be decomposed as 

(x, y) = \ {ip{x) + f/,(y)) a^h-Hx - y)) + a(x, y) (5.14) 



with a periodic function tp £ H 2 {C) satisfying 

\\iP\\ H i<C, U\\ H 2 < Ceh- 1 (5.15) 

and with 

Ikllffi < Ce- 2 h 3 . (5.16) 
More precisely, one has a = o~\ + 02 with 

\Wxfm ^ Ch ? ( 5 - 17 ) 

and with o~i of the form 

Mw) = \ + $(V)) oP{h-\x - y)) , 

where the Fourier transform of ip supported in {\p\ > eh^ 1 }. The Fourier transform 
of tp is supported in {\p\ < eh^ 1 }. 

We recall that the H 1 norm of an operator was introduced in (|3.7p . 

Proof. Step 1. We can write (|5.6p as 

0< [ (a(-,y)\K°(-ihV)-ah5(- ~ y)\a(- ,y)) dy < Ch 3 . (5.18) 



Here, the operator K^,(—ihV) acts on the x variable of a(x,y), and ( ■ | • ) denotes 
the standard inner product on L 2 (R). 

Now we recall that the operator — ahS(- — y) on i 2 (R) has a unique ground 
state, proportional to a°(/i _1 (- — y)), with ground state energy zero. There are no 
further eigenvalues and the bottom of its essential spectrum is A / tanh [^] > 2T. 
In particular, there is a lower bound, independent of h, on the gap. We write 



Mv)=(hj \a°(x)\ 2 dx) f a°{h- 1 {x~y))a{x,y)dx 



(5.19) 



and decompose 

a(x,y) = tp (y)a°(h^ 1 (x - y)) + cr (x,y) . 

Then (|5 . 18[) together with the uniform lower bound on the gap of Kj, — ahd( ■ — y) 
yields the bound ||oo||| < Ch 3 . We can also symmetrize and write 

a!(x,y) = a (x,y) + \ {^(x) - ^{y)) a\h~ l (x - y)) . (5.20) 

Then 

a{x, y) = \ (V>o(s) + Mv)) ^(h^ix - y)) + a x (x, y) (5.21) 

again with 

|| || I < Ch 3 . (5.22) 

This proves the first half of (|5.17|) . Before proving the second half in Step 4 below 
we need to study ip. 
Step 2. We claim that 

\M x )\ 2dx < C (5.23) 
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and 

/ \tp' {x)\ 2 dx < C. (5.24) 
Jc 

The first inequality follows by Schwarz's inequality 

II „.I|2 

2 



\Mx)\ 2 dx< 



hf R \cfl(x)\*dx 



and our bounds (|5.5|) and (|5.2p . In order to prove (|5.24[) we use again Schwarz's 
inequality, 

Lemma [2] below bounds the numerator by a constant times 

h- 2 f (a(-,y)\K°(-ihV)-ah5(- - y)\a{ ■ , y)) dy , 



and therefore (|5.24|) is a consequence of (|5.18p and f|5 . 2[) . 

Step 3. Next, we establish the remaining bound || Vcri ||| < Ch in (15.17|) . We 
use formula (|5.20|) for o\. First of all, using the fact that K{—ihV) > c(l - h 2 V 2 ) 
one easily deduces from (|5.18l) that HVooHl < Ch. Moreover, because of (|5 . 2|) and 
(EH 



/ W,{x)\ 2 I \a°(h-\x-y)){- 

JC JR 



dx dy < Ch . 



Finally, 



\Mx) - MvtfWih-H* ~ V))? dxdy 



Cxi 



= h~H \Mp)\ 2 f \(a )'{x)\ 2 sin 2 (}hpx) dx 

<hJ2 \P\ 2 \MP)\ 2 I \ (a Q )'{x)\ 2 x 2 dx<Ch, 

where we used (|5.24[) and the fact that J R |(a°)'(a;)| 2 x 2 dx is finite. This is a simple 
consequence of the fact that the Fourier transform of a is given by the smooth 
function 2 ^y/ 2 (^t) ■ This completes the proof of (|5.17[) . 

Step 4- Finally, for each e > h we decompose tJjq = ip + ip, where the Fourier 
transforms of ip and ip are supported in {|p| < eh^ 1 } and {\p\ > e/i -1 }, respectively. 
Clearly, the bounds ([533)1 and ((5341) imply ([535]) . 

Moreover, \\ip\\2 < Ce~~ 1 h and || 1 1 2 < C, and hence 

*2(x,v) = \ (#(x) + 4>{y)) a°(h'\x ~ y)) 
satisfies ||<72||fj-i < h 3 e~ 2 . This completes the proof of the proposition. □ 
In the previous proof we made use of the following 



1(5 



R.L. FRANK, C. HAINZL, R. SEIRINGER, AND J. P. SOLOVEJ 



Lemma 2. For some constant C > 0, 



xC 



\(V X + V y ) a(x, y)\ 2 dxdy 

<C j (a(- ,y)\K%(-ihV) - ahS(- - y)\a{- ,y)) dy (5.26) 

for all periodic and symmetric a (i.e., a{x,y) = a(y,x)). 
Proof. By expanding a(x,y) in a Fourier series 

a{x,y)= ]T e^+yy 2 a p (x-y) (5.27) 

and using that a p (x) — a p (—x) for all p £ 27rZ we see that (|5.26l) is equivalent to 

2 

K^(-ihV + hp/2) + K^{~ihV - hp/2) - 2a S(x/h) > -^h 2 p 2 (5.28) 

for all p 6 27rZ. This inequality holds for all p £ R, in fact, for an appropriate 
choice of C > 0, as we shall now show. 

Since Kj, > const(l + ft, 2 (— iV+p/2) 2 ), it suffices to consider the case of hp small. 
If k = Ao/ tanh [^1 > 2T denotes the gap in the spectrum of K^(-ihV) — ahd 
above zero, and h~ 1 / 2 <f>o(x/h) its normalized ground state, proportional to a°(x/h), 

K^(-ihV + hp/2) + K^(-ihV - hp/2) - 2aU 



> K 

> K 



e ihxp/2 (1 _ |0 O )^ Q |) e ~ihxp/2 + e ~ihxp/2 (1 _ |0 Q )(0 O |) e ihxp/2 



\Mx)\ 2 e- lhxp dx 



(5.29) 



In order to see the last inequality, simply rewrite the term as k(2 — |/)(/| — |<7)(s|), 
where |(/|g)| 2 = 1/ \4>o(x)\ 2 e~ lhxp dx\ , and compute the smallest eigenvalue of the 
corresponding 2x2 matrix. Since <^>o is reflection symmetric, normalized and sat- 
isfies J x 2 1 c/>o | 2 dx < oo (see Step 4 in the proof of Proposition [2]) , we have 



\Mx)\ 2 e~ lhxp dx 



This completes the proof of (|5.28l) . 



t>o(x)\ 2 (1 - cos(hpx)) dx > ch 2 p 2 



□ 



5.4. The lower bound. Pick a T with T BCS (T) < J" BCS (r°) and let ip be as in 
Proposition^ (depending on some parameter e > h to be chosen later). As before we 
let A(x) = -ip(x)A and define H A by (pTTj) . We also put T A = (1 +cxp( 1 SiJ A ))~ 1 - 
Our starting point is the representation 

^Bcs (r) _ j-BCS (ro) = -| Tr [ln(l + e"^) - ln(l + e^ Ho )] 
T 



+ — H(T,T A ) + A Re / ip(x)a(x, x)dx - ha I \a(x, x)\ 2 dx . (5.30) 
2 Jc Jc 

(Compare with (|5.7p .) According to the decomposition (|5.14|) which, in view of the 
BCS gap equation ()1.7|) . reads on the diagonal 



y(x, x) = ip(x)a°(0) + a(x, x) = — - + o~(x, x) , 

2ah 
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we can obtain the lower bound 

JTBCS (r) _ J-BCS (ro) >_ZTr [l n(1 + C -PH A) _ ln(1 + e -pH ^ 



T-L(T,T A ) ~ ah / \a(x, x)\ 2 dx . 
Jc 



For the first two terms on the right side we apply the semiclassics from Theorem [5] 
Arguing as in the proof of the upper bound and taking into account the bounds on 
ip from Proposition [5] we obtain 

^ BCS(r) _ J-BCS (ro) > ft 3 ( f GL W _ ^ _ ^3 



+ f H{T,r A )-ahJ \a{x,x)\ 2 dx . (5.31) 



Our final task is to bound the last two terms from below. In the remainder of this 
subsection we shall show that 

~W(r,r A )-a/i / \a(x,x)\ 2 dx > -C (eh 3 + e^h 4 ) . (5.32) 
Jc 

The choice e = h 1 / 3 will then lead to 

^BCS (r) _ J-BCS (ro) > ^3 ( £ GL _ ^ _ ^3+1/3 ^ 

which is the claimed lower bound. 

In order to prove (15.32)) we again use the lower bound on the relative entropy 
from Lemma [1] to estimate 



TH(T,T A ) > Tr 



^(r-TA) 2 



(5.33) 



tanh 2t^a 

The next lemma will allow us to replace the operator H A in this bound by Hq. 
Lemma 3. There is a constant c > such that for all sufficiently small h > 

" A >(l-ch)K^{-ihV)®l C 2. (5.34) 



tanh 

Proof. An application of Schwarz's inequality yields that for every < rj < 1 

Hi > (1 - n) (H A f - ^(AlU - 1WI + h 2 \\W\\l) . 
The expansion formula |14( (4.3.91)] 

2 + E( 2 - 



tanh(x/2) x 2 /4 + k 2 TT 



shows that x i— > \fxj tanh yfx is an operator monotone function. This operator 
monotonicity implies that 



, (1 - v)- 1/2 VH'i + V-'i^U ~ l|| 2 oc + h 2 \\WPJ 



Kj,(-iKV) ® Ic 2 < 



tanh -L(l - rj)-V*y/E* + V^iH - MIL + ^WK) 



2T ^ 

<(1 _ r i/ 2 y/Hj + y-^W-nio + hVlll,) 

tanh^^A + ^A^P - 1||L + ^ 2 II^||L) 



<(i-^(i + ^(A8^-iBi + AW.)) 5 ^ 
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for < 77 < 1. The Sobolev inequality and (|5.15j) show that |] VII 00 < C||i/>||ffi < C, 
and hence the lemma follows by choosing 77 = h. □ 

To proceed, we denote a a = [Ta] 12 and recall from Theorem [3] that 
«a = j (ipK^-ihV)- 1 + K^i-ihV)- 1 ^) +771+772 

= 2 (^ a ° + a< V) +V1+V2 

with 771 and 772 satisfying the bounds (|3.9[) and (|3 . 10|) . The second equality follows 
from the explicit form (|1.10[) of a°. Comparing this with (|5.14l) we infer that 

a = + cr — f]i — 1J2 ■ (5.35) 

Then (^33j) and ([534]) imply that 

■j H(T, Ta) — ah / \a(x, x)\ 2 dx 
Jc 

> (1 - ch) Tr K%(-ihV)(a - a a) (a - oa) - ah J \a(x, x)\ 2 dx 

> (1- ch)TrK%(-ihV)<ra- ah J \a(x, x)\ 2 dx 

- (1 - ch)2ReTr K%(-ihV)a(rjr+rfr) ■ (5.36) 

In order to bound the first term on the right side from below we are going to 
choose a parameter p > such that ch + p < 1/2. Here c is the constant from 
(15.341) . (Eventually, we will pick either p = or p = 1/4, say.) Note that 

(1 - ch)TrK%,(-ihV) - ahS = pK^(-ihV) + (1 - 2ch-2p)(K%,(-ihV) - ahS) 

+ (ch + p)(K^(~ihV) ~ 2ahS) . 

We recall that the operator K^(—ihV) — ah5 is non-negative and that the operator 
K^(—ihV) — 2ahd has a negative eigenvalue of order one (by the form boundedness 
of 5 with respect to K^,(—iV)). Hence K^(—ihV) — 2ah5 > —C\ with a constant 
C\ independent of h. (In the following it will be somewhat important to keep track 
of various constants, therefore we introduce here a numbering.) Moreover, using 
the fact that K%(-ihSJ) > ci(l - /i 2 V 2 ) we arrive at the lower bound 

(1 - ch) Tr K%(-ihV) - ahS > c x p{\ - h 2 V 2 ) - d(ch + p) , 

which means for the first term on the right side of (I5.36P that 

(l-ch)TrK%(-ihV)aW-ah J \<r(x, x)\ 2 dx > cxp\\a\\ 2 H i — d(cft.+/o)||cr||| . (5.37) 

We now turn to the second term on the right side of (|5.36p . Theorem [31 together 
with the bounds (|5 . 1 5 1) on ip, implies that \\r]i + ^W^i < Ce 2 h 3 . This bound, 
combined with ||ct||^-i < Ce~ 2 h 3 from Lemma [U however, is not good enough. (It 
leads to an error of order h 3 .) Instead, we shall make use of the observation that 
in the decompositions a = oi + <7 2 and 771 + 772 one has 

Tri^(-iW)o- 2 ^r = 0. 

This can be seen by writing out the trace in momentum space and recalling that 
the Fourier transform of the ip involved in 01 has support in > e/i -1 }, whereas 
the one of the i\> involved in 771 has support in {|p| < eh^ 1 } (see also (|6.22[1 ). 
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Using the estimates (|5.17|) and (|3.10j) on u\ and r\ 2 we conclude that 
|Trjr£(-iW)cr(?7i +772)] < |Tr K%{-ih V)<7i)yT| + |Tr Jf§,(-iW)cr^| 

< C 2 (eh 3 + h 5 ' 2 \\a\\ H ^j . (5.38) 
Combining ([536]) . (1537)) and (f5T38|) we find that 
-j "H(r, Ta) — a/i y |cr(a;, x)\ 2 dx 

> dPlWWm - Ci{ch + p)\\a\\l - 2C 2 (eh 3 + h 5 / 2 \\a\\ H ^j . (5.39) 

Next, we are going to distinguish two cases, according to whether 4Ci||(t||| < 
ci||cr||^i or not. In the first case, we choose p = 1/4 and h so small that ch+p < 1/2. 
In this way we can bound the previous expression from below by 

idlHI^ - 2C 2 (eh 3 + h^ 2 \\a\\ H i^ > Sc^Cih 5 - 2C 2 eh 3 . 

This proves the claimed (indeed, a better) bound (|5.32[) in this case. 

Now assume, conversely, that 4Ci ||cr|| § > ci||er||^i. Then we choose p = and 
bound (|5.39|) from below by 

- C lC h\\a\\l - 2C 2 (eh 3 + h 5 ' 2 \\a\\ H ^j 

> -C lC h\\a\\l - 2C 2 (eh 3 + h b ' 2 (4Ci/ Cl ) 1/2 ||a|| 2 ) . 

The bound (|5 . 16[) on ||er|| 2 now leads again to the claimed lower bound (|5 -32[) . 

This concludes the proof of the lower bound to the free energy in Theorem [T] 
Concerning the statement about approximate minimizers we note that J 7 (F°) — 
J- B (To) = 0(h 3 ) and that our a-priori bounds on a in Proposition [2] remain true 
under the weaker condition that J rBCS (r) < J" BCS (r°) + Ch 3 . We leave the details 
to the reader. 

6. Proof of semiclassical asymptotics 

In this section we shall sketch the proofs of Theorems [5] and [3] containing the 
semiclassical asymptotics. We shall skip some technical details and refer to [7] for 
a thorough discussion. 

6.1. Preliminaries. It will be convenient to use the following abbreviations 

k = -h 2 V 2 - v + h 2 W(x), k a = -h 2 V 2 -p. (6.1) 

We will frequently have to bound various norms of the resolvents (z — fc) _1 for z 
in the contour T defined by Imz = ±7r/(2/3) for j3 > 0. We state these auxiliary 
bounds separately. 

For p > 1, we define the p-norm of a periodic operator A by 

\\A\\ p =(Tv\A\n 1/p (6.2) 

where Tr stands again for the trace per unit volume. We note that for a Fourier 
multiplier A{— ihV), these norms are given as 

\\A(-ihV)\\ p = h-^ \A(q)f^j . (6.3) 
The usual operator norm will be denoted by ||^4||oo- 
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Lemma 4. For z — t ± i7r/(2/3) and all sufficiently small h one has 

\< u\-i\\ s-nu-lh f ^ 1/(2p) /ort»l .... . .„ 

\{z-k) l \\ n <Ch VP x | |t| _ 1+1/(2p) } ^ << _ 1 */l<P<oo, (6-4) 



as well as 

Ik^HL^xj^ . (6.5, 

Proof. The estimates are easily derived with fco instead of k by evaluating the 
corresponding integral. Since the spectra of k and ko agree up to 0{h 2 ) the same 
bounds hold for k. □ 

6.2. Proof of Theorem[2l The function / in (|3.2p is analytic in the strip | Imz| < 
7r, and we can write 



f{pH A )-f{[3H ) = ^-.Jf{f3z 



1 1 



z - H A z — H Q 



dz , 



where T is the contour z — r ± i-^, r 6 R. We emphasize that this contour 
representation is not true for the operators /(/3i?A) and f(/3Ho) separately (because 
of a contribution from infinity), but only for their difference. 
We claim that 



[fWH A )} n - [f(f3H A )} 22 - [3[H A ] 22 . (6.6) 

Recall that [ • ]ij denotes the ij element of an operator- valued 2x2 matrix. To see 
(I6.6[) . we introduce the unitary matrix 



U = 

and note that 



1 

-1 



[f(0H A )] n = - [Uf(f3H A )U} 22 . 
On the other hand, UH A U = —H A , which implies that 



Uf{PH A )U = f(-j3H A ) = f(-[3H A ) . 

The claim (|6.6[) now follows from the fact that f(—z) = f(z) — z. 

Subtracting (I6.6[) and the corresponding formula for H and noting that H A and 
H coincide on the diagonal we find that the two diagonal entries of f((3H A ) — 
f(/3H ) are complex conjugates of each other. Since their trace is real we conclude 
that 



Tr [f(pH A ) - f(f3H )} = - f f((3z) Tr 



1 1 



z - H A z- H 



dz . 

li 



(For technical details concerning the interchange of the trace and the integral we 
refer to [7].) 

The resolvent identity and the fact that 
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is off-diagonal (as an operator-valued 2x2 matrix) implies that 

2" 



Tr 



1 



1 



z - Ha z - H 
1 



li 



Tr 



Tr 



1 



1 



1 



z - H \ z- H 

4 



z — Ho V z — Hi 



Ti- 



ll 



1 



z — H a \ z — H 



J n 



=:h+I 2 + I 3 . 
In the following we shall prove that 



1 

7TZ Jp 



f(/3z)h dz 



hp 2 



dq 

2^ 



+ ^rVlll / ( 5 i(/3(g 2 -M)) + 2/3g 2 .g 2 (/3(g 2 -M))) 



+ / gi(P(q 2 -ri) 



dq 
27 



o(h 5 )U\\ 2 H2 



dq 
2^ 



(6.8) 



- / f{fiz)h 



dz = 



h 3 f3 3 



mi 



gi(/3(g 2 - //)) dq 
q 2 — ji 27r 



+ 0(h 5 )W\\ 3 H1 W\ H 2 (6.9) 



and 



1 



f(Pz)I 3 dz = 0(h 5 )M 6 Hl 



(6.10) 



This will clearly prove (|3.6p . We will treat the three terms 1$, I2 and 1\ (in this 
order) separately. 

I3 : With the notation k introduced in (I6.ip at the beginning of this section, we 
have 



h = Tr 



1 



z-H A 



A — 3— A t — 3— A — ^— A t — 5— A — 5— A t — — 



11 z + fc z — k z + k z — k z + k z — k 

Using Holder's inequality for the trace per unit volume (see [7]) and the fact that 
\z - H A \ > n/(2/3), we get 

i/3i<^ 8 w2oii(*-fc)- i n;ii^+*)- i K. 

Together with (|6.4[) . this yields 

Here it was important to get a decay faster than |z| -2 , since we need to integrate 
^3 against the function / which behaves linearly at —00. Since Halloo < C||^>||jyi 
by Sobolev inequalities we have completed the proof of (|6.10[) . 



I2 : We continue with 



h = Tr — ^— -A — ^— -A^ — 5— A — ^— rA^ — — 



z — k z + k z — k z + k z — k 
By the resolvent identity we have 

1 1 



z — k z — ko z — k { 



' h 2 W— 



(6.11) 
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Using Holder as above, we can bound 
1 1 



Tr 



fe 



A 1 At 1 A 1 At 1 

A -A 1 -A -A'- 



3.12) 



z + k z — k z + k z — k 

By IJB3D and this is bounded by Ch 5 \\ip\\%i, c) (l + |*| 5/2 ) _1 . What we effec- 
tively have achieved for this error is, therefore, to replace one factor of (z — fc) -1 in 
I 2 by a factor of (z — fco) -1 

In exactly the same way we proceed with the remaining factors (z — fc) _1 and 
(z + fc) _1 in The only difference is that k might now be replaced by kg in the 
terms we have already treated, but this does not effect the bounds. 

The final result is that (ni) -1 Jr fifi z ) h dz equals 



1 



JWTr 



1 



1 



-A+- 



1 A 1 At 

-A — A*- 



1 



z — fco z + ko z — ko z + fco z — fco 

and it remains to compute the asymptotics of the integral. 

Let us indicate how to perform the trace per unit volume Tr[ 
integrals the trace can be written as 



dz + 0{h 5 )\\ 



Ih 1 i 



In terms of 



(27T) 



dx\ I dx2 

e ipi(x 1 -x 2 ) 



dx 3 I dx 



dp 1 dp 2 dp 3 dp i i/)(xi)ip(x2)ip(xa)ip(x4) 



IP2(X2-X 3 ) 



oiP3(x 3 -x 4 ) 



ipi{xi-xx) 



(6.13) 



(z - (h 2 p\ - fij) 2 z + {h 2 p\ - n) z— (h 2 pl - n) z + (h 2 pl - fi) 
Since ip is periodic with period one we have 

E HhY x > h - 

We insert this into the above integral and perform the integrals over x 2 , x 3 , X4. This 
leads to ^-distributions such that we can subsequently perform the integrals over 
P2,P3,P4, as well as the integral over x\. In this way we obtain 



1 



1 



1 



-At- 



1 



-A- 



1 



-At- 



1 



dz 



^ E 

Pi,P2,P3£2ttZ 



z — fco z + fco z — ko z + fco z — fco 

ip(pi)^*(P2)^(P3)^*(-pi -P2 - p 3 )F{hpi,hp 2 ,hp 3 ) 



with 



F(Pl,P2,P3) 



X 



dzfWz) 



dq 



1 



1 



2tt (z - (3((q + Pl + p 2 + p 3 ) 2 + fj)Y 

1 1 



z + f3{{q+ Pl +p 2 ) 2 -fi) z - /3((q + pi) 2 - n) z + f3(q 2 -v) 
The leading behavior is given by 

F(0,0,0) E kpi)r(P2Mp 3 )r(-Pi-P2-P3) = F(Q,Q,0)U\\j. 

pi,P2,P3&2irZ 

The integral -F^O, 0, 0) can be calculated explicitly and we obtain 

/? 3 f 9lW(q 2 -^)) dq 



F(0,0,0) 



2tt 
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with gi from (|3.4j) . In order to estimate the remainder we use the fact that |7] 

|^(P1,P2,P3)--F(0, 0,0)| < const (pl+pl+p 2 3 ) . 
Using Schwarz and Holder we can bound 

Pi fripiWfawipaWi-Pi ~P2 -pa) < const H^II^II^IIh! 

and equally with p\ replaced by p\ and p\. Hence we conclude that 

— f(/3z)I 2 dz = h 3 V ${px)^(p2)$(pz)&{-pi-P2-pz)F{hpx,hp2,hpz) 
m J r * — ' 

pi,P2,P3£27rZ 

+ 0(^)11^11^ 
= A 3 f (0,0,0)11^111 + 0(h 5 )Mf Hl \\i>\\ H 2 . 
This is what we claimed in (16.91). 



Ii : Finally, we examine the contribution of 
h = Tr 



1 A 1 At 1 

A -A'- 



z — k z + k z — k 
Using the resolvent identity (|6.11[) we can write Ii = 7" + l\ + If, where 

I{ 1 = Tr 



1 A 1 At 1 

-A — AT- 



and 



II = Tr 



z — ko z + fco z — ko 



' (fc-fc )— ^A— ^A*- ' 



z- k 



z — ko z + ko z — ko 



' -A — ^—{k - k) — i-^A t - ' 



z — ko z + fcg 



z + fco z — ko 



z — ko z + fco z — ko z — fco 



The part If consists of the rest. We claim that 



1 

ni 



f{f}z)iidz = -^w\\i 1 <,n(.;(,r'-/<n 



dq 
2^ 



U'Wl / (gi({3(q 2 -n))+2j3q 2 g 2 {P( q 2 - fJ ,))) 



dq 
2^ 



0(h 5 )\\ 



\H 2 



3.14) 



m 
and 



h 3 p 3 



- I f(f3z)I b 1 dz=^{4>\W\il>) I gi (f3( q ^^)^ + 0(h b )yj\\ H ,U\\ m (6.15) 



dq 



2tt 



- / f^z)ifdz=o(h")ur H1 . 



(6.16) 



Clearly, this will imply (|6.8I) . 

We begin with If. These terms contain at least five resolvents, where at least 
two terms are of the form (z — fc#) _1 and at least one term of the form [z + fc#) _1 . 
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(Here k# stands for any of the operators k or ko.) Moreover, they contain at least 
two factors of k — ko. The terms are either of the type 

A = Tr^—(k-k )^—A— L-iko-k)^—^^— (6.17) 

z — ko z — k z + ko z + k z — ko 

(at least three minus signs) or of the type 

B = Tr — i— (fc — ^O - 7-77(^0 — TT^ — ~~T ( 6 - 18 ) 

z — ko z + ko z + ko z + k z — k 

(only two minus signs) . Terms of the first type we bound by 

|A|<c/ l 6 ||^||LllV-llLll(^-fco)- 1 llLll(^ + M-Nl3ll(^ + fc)- 1 || 3 ||(z-fc)- 1 ||3. 

By (RO|) and (JOJ) this can be estimated by C7i 5 |z|~ 2+1/6 if Rez > 1 and by 
C/i 5 |z| -2-1 / 6 if Rez < —1. This bound is finite when integrated against f((3z). 

Terms of type B can be bounded similarly by replacing z by —z. Indeed, we 
note that since J r z B dz = 0, we can replace f(/3z) by f(—/3z) — ,f{(3z) — fiz in 
the integrand without changing the value of the integral. I.e., we can integrate B 
against a function that decays exponentially for negative t and increases linearly 
for positive t, instead of the other way around. These considerations lead to the 
estimate (|6.16[) . 

Next, we discuss the term After doing the contour integral the term gives 
(m)- 1 f }{Pz)I a l dz = h H(p)\ 2 G(hp) 

with 

P f tanh (\P({q + p) 2 - /i)) + tanh (±f3(q 2 ~ fi)) dq 



aip)-U 

L JR 



[p + q) 2 + q 2 -2fi 2?r 
By definition (|3.3p we have 



G(0) = -^ / goifitf-Mp. 

2 Jm. 2lT 



Integrating by parts we can write 

G"(0) = ^ jf ( 9l (f3(q 2 M)) + m 2 92(P(q 2 - M))) ^ 

with (71 and gi from (|3 .4[) and (|3.5p . Moreover, one can show that [7] 

\G(p) - G(0) - \p 2 G"{0)\ < Cp 4 . 
From this we conclude that 

i Jj{Pz)I a 1 dz = h Y \m\ 2 (G(0) + ±G"(0)h 2 p 2 )+O(h 5 )U\\ 



2 



hG{<S)U\\l + hG"(o)h 3 W\\i + o(h 5 )M 2 



p£2-nZ 

■"- 1 \V)">~ WW II 2 t- )wip\\ fJ , 
which is what we claimed in (|6.14|1 . 

Finally, we proceed to l\. After the contour integration we find 



\ [ f{pz)I h 1 dz = h? r(pM<l)W(-p-q)L(h P ,hq), 
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where 



L(p,q) = f3 3 / L(p,q,k) 



dk 
2^ 



with 



L{p,q,k)=-. f In (2 



+ 



1 



1 



1 



z + k 2 ~(iz — p 2 +[iz — q 2 +p 
1 1 



z — p 2 + /! z — q 2 + fi z + k 2 — fi 



dz . 



We have 



£(o,o) 



P 3 



giW(k 2 -p)) 



dk 
2^ 



and (see [7] for details) 

\L(p,q)-L(0,0)\<C(p 2 +q 2 ) . 
By the Schwarz inequality we can bound 

£ \r( P H(q)W(-p-q)(p 2 + q 2 )\ < C\\W\\ 2 W\\h*W\\b* , 

and obtain 

- [ f((3z)I b 1 dz = h 3 L(0 7 0) V ^(-pWg)W(p-g) + 0(^)11^11^11^11^ 



/j 3 /3 3 



dk 



WW]*) / S i(/3(A i -M))^ + 0(/i 5 )||V||^|IV'll^ 



2tt 



This concludes the proof of Theorem [51 



6.3. Proof of Theorem [3l Since the function p in (|3.8|) is analytic in the strip 
| Iva.z\ < 7r, we can write [p{j3H/^)] 12 with the aid of a contour integral representa- 
tion as 

1 



WH A )} 12 = f p(fiz) 
2m J r 



z-H A 



dz , 



(6.19) 



12 



where T is again the contour Imz = ±7r/(2/3). We expand (z — Ha) 1 using the 
resolvent identity and note that, since Ha = Ho + S with a Hq diagonal and 5 off- 
diagonal, only the terms containing an odd number of S's contribute to the 12-entry 
of (z — Ha) ■ In this way arrive at the decomposition 



p(/3H A )] 12 = Vo + V1+V2+V21 



where 



h 

Am' 



p{Pz) V 



+ 72— T2 ^ I dz ■ 



m = hj v p ^ z) {j^v ^ ko] ^h 2 + ^ fco] ir^) dz ' (6 



'/2 



27TJ 



p(/3z) 



2T- fe 



W -tp + ip — — w 

z — k z + fcn 



2 + fc 



(6.20) 

(6.21) 
22) 
(6.23) 



and 



4 



'2m 



1 i^jh^jh^^ 



z-H A 



dz . 



(6.24) 
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A simple residue computation yields 



"4^ 



p(0k o ) - p(-0ko) , ptfko) - p(-0ko) 



k 



1> 



hd 



(ipg (Pk ) + g (/3k ) V) , 



which is the main term claimed in the theorem. In the following we shall prove 
that 

\\mfm<Ch 5 Uf H2 , 
W\ 2 m<Ch 5 U\\ 2 m , 

and 

\\v b 2\\ 2 m<Ch 5 \\it m . 
This clearly implies Theorem [3] 



(6.25) 
(6.26) 

(6.27) 



7/1 : The square of the H 1 norm of 771 is given by 

\\m\\h=h J2 \Mp)\ 2 J(h P ) 

with 

?4 



J(P) = ^T [ ((q + P) 2 - q 2 ) 2 (l 



V) \F(q+ P ,q)-F(q,q 



and 
F(p,q) 



p 2 - /i 1 + ef(P 2 -^ 1 + 



n ,2 d<? 
P)l 27 



= /3(p 2 + 9 2 -2 M ) 



1 



q 2 - 2p 



One can show that < J(p) < Cp 4 [7], which yields the desired bound ()6.25j) . 



r/2 : This term is a sum of two terms and we begin by bounding the first one, 
that is, ~h 3 (2Tri)- 1 j p((3z)(z - ko^Wiz - fc)" 1 ^ + k)- 1 dz. Using Holder's 
inequality for the trace per unit volume we find that the square of the H 1 norm of 
the integrand can be bounded by 



Tr 



1 - h 2 V< 



w 



1 



^7 



1 



\z - k \ 2 z-k \z + k\ 



- k 



W 



< 



1 - h 2 V 2 



- 1 1 1 2 
2 • 



In order to bound this we use (|6.4I) and (|6.5p . as well as the fact that ||(1 — h 2 V 2 )\z — 
k \- 2 \\oo is bounded by C|z| _1 if Rez < -1 and by C\z\ if Rez > 1. (This follows 
similarly as ()6.5[) .) In particular, we conclude that for Rez < — 1 the previous 
quantity is bounded by C/i -1 !!?/;!!^^! -7 / 2 . The square root of this is integrable 
against p((3z) and we arrive at the bound Ch b / 2 HV'lloo for the H 1 norm. For the 
positive z direction, we notice that p((3z) decays exponentially leading to a finite 
result after z integration. 

For the second term in vfe we proceed similarly. It is important to first notice 
that p(z) = 1 — p{— z), however, and that the 1 does not contribute anything but 
integrates to zero. Proceeding as above we arrive at (|6.26[) . 
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772 : Finally, we consider n^. Using Holder's inequality for the trace per unit 
volume and bounding [(z — Ha) ^22 by 2/3/ir for z € T we find that the square of 
the H 1 norm of the integrand is bounded by 



uz-krmuz+krm 

00 

Similarly as in the bound for ?yf one can show that for Rez < — 1 this is bounded 
by ChT^Wl^z]- 7 ' 2 . This leads to ^?Fj) . 
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